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JEAN BERARD 

Abstract. We consider a branching-selection particle system on TL with A*' > 1 
particles. During a branching step, each particle is replaced by two new particles, 
whose positions are shifted from that of the original particle by independently 
performing two random walk steps according to the distribution p5i + (1 — p)5o, 
from the location of the original particle. During the selection step that follows, 
only the A rightmost particles are kept among the 2A particles obtained at the 
branching step, to form a new population of N particles. After a large number 
of iterated branching-selection steps, the displacement of the whole population of 
N particles is ballistic, with deterministic asymptotic speed vn{p)- As N goes 
to infinity, vn{p) converges to a finite limit Voo{p)- The main result is that, for 
every < p < 1/2, as N goes to infinity, the order of magnitude of the difference 
Vcx,{p) — ujv(p) is log(A)~^. This is called Brunet-Derrida behavior in reference to 
the 1997 paper by E. Brunei and B. Derrida "Shift in the velocity of a front due to 
a cutoff' (see the reference within the paper), where such a behavior is established 
for a similar branching-selection particle system, using both numerical simulations 
and heuristic arguments. The case where 1/2 < p < 1 is briefly discussed. 



1. Introduction 

In [ll[5], E. Brunet and B. Derrida studied a branching-selection particle system 
on Z enjoying the following property: as the number N of particles in the system 
goes to infinity, the asymptotic speed of the population of particles in the system 
converges to its limiting value at a surprisingly slow rate, of order log(A^)~^. This 
behavior was established both by direct numerical simulation of the particle system, 
and by mathematically non-rigorous arguments of the following type: having a finite 
population of particles instead of an infinite number of particles should be more or 
less equivalent, as far as the asymptotic speed is concerned, to introducing a cutoff 
value of e = in the deterministic equations that govern the time-evolution of 

the distribution of particles in the infinite-population limit. In turn, these equations 
can be viewed as discrete versions of the well-known F-KPP equations, and the 
initial problem is thus related to that of assessing the effect of introducing a small 
cutoff in F-KPP equations, upon the speed of the travelling wave solutions. In 
turn, this problem was studied by heuristic arguments and computer simulations 
(see in [5])) and rigorous mathematical results for this last problem have recently 
been obtained, cite OIUIH]. Another approach is based on adding a small white- 
noise term with scale parameter o"^ = in the Fisher-KPP equation, see [6], and 
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rigorous results have been derived for this model too, see [3 EH]- However, to our 
knowledge, no rigorous results dealing directly with a branching-selection particle 
system such as the one originally studied by Brunet and Derrida, are available. In 
this paper, we consider a branching-selection particle system that is similar (but 
not exactly identical) to the one considered by Brunet and Derrida in [H [5] . To be 
specific, we consider a particle system with N particles on Z defined through the 
repeated application of branching and selection steps defined as follows: 

• Branching: each of the A'^ particles is replaced by two new particles, whose 
positions are shifted from that of the original particle by independently per- 
forming two random walk steps according to the distribution p6i + {1 — p)6q, 
from the location of the original particle; 

• Selection: only the N rightmost particles are kept among the 2N obtained 
at the branching step, to form the new population of A'^ particles. 

In Section [3l it is proved that, after a large number of iterated branching-selection 
steps, the displacement of the whole population of N particles is ballistic, with deter- 
ministic asymptotic speed vn{p), and that, as N goes to infinity, vn{p) increases to a 
finite limit Voo {p) (which admits a more explicit characterization) . The main results 
concerning the branching-selection particle system are contained in the following 
two theorems: 

Theorem 1. For every < p < 1/2, there exists < C*(p) < +cxd such that, for 
all large N, 

(1) Vooip) - vn{p) > CM log(iV)"' 

Theorem 2. For every < p < 1/2, there exists < C*{p) < +oo such that, for 
all large N , 

(2) Vooip) - vn{p) < C*{p) log(Af)-2 

In the case p = 1/2 or 1/2 < p < 1, the behavior turns out to be quite different, 
as stated in the following theorems. 

Theorem 3. For p = 1/2, there exists < C*(l/2) < C*(l/2) < +oo such that, 
for all large N , 

(3) C,{1/2)N-' < vM - Mp) < C*{l/2)N-\ 

Theorem 4. For p > 1/2, there exists < C*(p) < C*{p) < +oo such that, for all 
large N , 

(4) C,{p)N < - logMp) - vn{p)) < C*{p)N. 

The rest of the paper is organized as follows. In Section O we provide the precise 
notations and definitions that are needed in the sequel. Section [3] contains a discus- 
sion of various elementary properties of the model we consider. Section [4] contains 
the proof of Theorem [H while Section [5] contains the proof of Theorem [21 The proofs 
of Theorems [3] and H] are sketched in Section [H Section [7] contains some concluding 
remarks. 
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2. Notations and definitions 

Throughout the paper, p denotes a fixed parameter in ]0, l/2[. Since the particles 
we consider carry no other information than their position, it is convenient to repre- 
sent finite populations of particles by finite counting measures on Z. For all A'^ > 1, 
let Cjv denote the set of finite counting measures on Z with total mass equal to N , 
and C the set of all finite counting measures on Z. 

For G C, the total mass of v (i.e. the number of particles in the population 
it describes) is denoted by M(z^). We denote by max;/ and minz^ respectively the 
maximum and minimum of the (finite) support of jJL. We also define the diameter 
^(1^) := max v — min u. Given two positive measures v on Z, we use the notation 
/X < to denote the fact that ^(x) < i^(a;) for every x G Z. On the other hand, we 
use the notation -< to denote the stochastic order between positive measures: -< u 
if and only if ijl{[x, +co\) < ^'([x, +oo[) for all x G Z. In particular, ^ -< v implies 
that M(/x) < M{u), and it is easily seen that, if = and = Ylil^i^ ^yt^ 

with xi > • • • > xa^(^) and yi > • • • > yM(u)i -< is equivalent to Xj < yi for all 
1 < f < M{v). From the order -< on C, we define the corresponding stochastic order 
on probability measures on C and denote it by -<^: given two probability measures 
Q and i? on C, Q -« R means that for every bounded non-decreasing function 
/ : (C, -<) M, Q{f) < R{f)- An equivalent definition is that there exists a pair 
of random variables {X, Y), X and Y taking values in C, such that X Q,Y R, 
and X -<Y with probability one. 

In this context, the dynamics of our particle systems can be defined through the 
following probability kernels. Let us first define the branching kernel p^^' on Cat xC2Ar 
as follows. Given v = Yld=i^xi £ Cat, pff'{v,-) is the probability distribution of 

I]i=i ^xi+Yi^i + Sxi+Y2,i G C2N, where (l£,i)i<i<Af, e=i,2 is a family of i.i.d. Bernoulh 
random variables with common distribution p6i + (1 — p)So. Then, we define the 
selection kernel pff'" on C2N x Cn as follows. Starting from v = Ylf^i G C2N, 
where xi > ••• > X2n, P^^'iy-,-) is the Dirac distribution concentrated on the 
counting measure Yl!i=i ^xi ■ 

The kernel pN on Cn x Cn governing the evolution of particle systems with TV 
particles is then defined as the product kernel pN := pff'pff'"- 

In the sequel, we use the notation {X.^)n>o to denote a Markov chain on Cat 
whose transition probabilities are given hy pN, and which starts at Xq := N60. We 
assume this Markov chain is defined on a reference probability space denoted by 

Let ~jv denote the equivalence relation on Cjv defined by ~Ar /i if and only if 
there exists m G Z such that u is the image measure of jj, by the translation of Z 
X ^ X + m. Let ttn denote the canonical projection from Cn to Cn/ ""^n- 

3. Elementary properties of the model 
Proposition 1. For all I < Ni < N2, and n G Cni and v G Cn^ such that n < f, 
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Proof. The proof is by coupling. Let fi = Yli=i ™^ ^ — Z]j=i ^Vi ' with 
xi > ■ ■ ■ > xatj and y\ > ■ ■ ■ > such that Xi < in for all 1 < i < A''i. Let 
(^,i)i<j<Af2, £=1,2 denote a family of i.i.d. Bernoulli random variables with common 
distribution p6i + (1 — p)6q. By definition, the couting measure defined by ^Br. '■= 
Ei=i ^x,+Y^,,+Sx,+Y2,, has the distribution fipff[-, and usr. ■= Ei=i ^y,+Yi,^+Sy,+Y2,, 
has the distribution np^^^. It is easily checked that hbt. -< i^Br., owing to the fact that 
Xi + Ye^i <yi + Yi^i for all 1 < i < iVi and £ = 1, 2. We deduce that fipfi[' ~<~< f^P^l'- 
Then, it is obvious from the definition that pff_l' and pff^' preserve -<~<. □ 



Proposition 2. For all N >1, and all n > 0, d{X^) < with probability 



Proof. Let m := [ iog(2) "l"^'^' "^^^ result is obvious forn = 0, . . . , m, since we perform 
or 1 random walk steps starting from an initial condition where all particles are 
at the origin. Now consider n > m, and let y = maxX^^. Assume first that 
minX^ < y for alln + 1 — m < k < n. Since all the random walk steps that 
are performed during branching steps are > 0, this implies that all the particles 
descended by branching from a particle located at y at time n — m, are preserved 
by the selection steps performed from X^_^ to X^ . Since there are 2^" > A'" such 
particles, this is a contradiction. As a consequence, we know that there must be an 
index n + l — m<k<n such that min X^ > y. Since by construction t i— > min X[^ 
is non-decreasing, we deduce that min X^ > y. Now, since all the random walk 
steps that are performed add or 1 to the current position of a particle, we see 
from the definition of y that maxX^ < y + m. As a consequence, d{X^) = 
maxX;^^ — minX^^ <m. □ 

Proposition 3. For all N > 1, the kernel p^ is compatible with the canonical 
projection vrjv, that is: for all z^i,z^2 G Cn such that vi ~iv V2, o,nd all ^ G Cn/ ^n, 

Proof. Immediate: everything in the definition of the branching and selection steps 
is translation- invariant. □ 

Proposition 4. For all N > 1, the sequence TTN{Xn)n>o is an ergodic Markov 
chain on a finite subset of Cn/ ^n- 

Proof. The fact that '^N{^n)n>o forms a Markov chain on C^/ ~Ar is an immediate 
consequence of Proposition [3l Let Sn '■= G Cn/ ~Ar; 3n > 0, P(7rAr(X^) = 
^) > 0}. From Proposition [2l we see that Sn is in fact a finite set. On the other 
hand, given any ^ € Sn, it is quite easy to find a finite path in Sn starting at ^ and 
ending at iTNiXdo) that has positive probability. As a consequence, the restriction of 
TTN{Xn)n>o to Sn is an irreducible Markov chain. As for aperiodicity, the transition 
TTNiNSo) I— > ttn{X5o) has e.g. a positive probability. □ 

Corollary 1. There exists < Wiv(p) < +oo such that, with probability one, and in 



one. 
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lim n minX„ = lim n maxX 
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Proof. Note that, in view of Proposition O if the two hmits in the above statement 
exist, they must be equal. Then observe that, for all n > 0, conditionally upon X:^ , 
the distributions of the increments max — max and min -^^i — min X^ 

depend only on itj<i{X^). The result then follows by a classical argument using the 
law of large numbers for additive functionals of ergo die Markov chains. □ 

Proposition 5. The sequence {vn{'p))n>i is non- decreasing. 

Proof. Consequence of the fact that, when A''i < N2, N160 -< N2S0, and of the 
monotonicity property [H □ 

We can deduce from the above proposition that there exists < Voo{p) < +00 such 
that limjv^+oo vv(p) = iJoo{p)- A consequence of the proofs of Theorems [H and [2] 
below is that Voo(p) is in fact equal to the number v{p) characterized as the unique 
root of the equation A(x) = log(2), where x € [0, 1] is the unknown, and where A is 
the large deviations rate function associated with sums of i.i.d Bernoulli(p) random 
variables, i.e. A{x) := x log(x/p) + (1 — x) log(^^) for x G [0,1]. We note for future 
use that p < v{p) < 1 since p s]0, l/2[. 

4. The upper bound 

The essential arguments used here in the proof of the upper bound, are largely 
borrowed from the paper [TT] by R. Pemantle, which deals with the closely related 
question of obtaining complexity bounds for algorithms that seek near optimal paths 
in branching random walks. In fact, the proof of Theorem [T] given below is basically 
a rewriting of the proof of the lower complexity bound in [11] in the special case 
of algorithms that do not jump, with the slight difference that we are dealing with 

independent branching random walks being explored in parallel, rather than a 
single branching random walk. 

To explain the connexion between our model and the branching random walk, con- 
sider the following model. Let BRWi, . . . ,BRW7v denote independent branching 
random walks, each with value zero at the root, deterministic binary branching, and 
i.i.d. displacements with common distribution p6i + (1 — p)So along each edge. For 
1 < i < iV, and n > 0, let BRWi(n) denote the set of vertices of BRWj located at 
depth n in the tree, and let T„ := BRWi(n) U • • • U BRWAr(n). For every n, fix an 
a priori (i.e. depending only on the tree structure, not on the random walk values) 
total order on T„. We now define by induction a sequence (Gn)n>o of subsets such 
that, for each n > 0, Gn is a random subset of with exactly N elements. First, 
let us set Go := Tq. Then, given n > and G„, let Hn denote the subset of T{n + 1) 
formed by the children of the vertices in Gn- Then, define Gn+i as the subset of Hn 
formed by the vertices that are associated with the largest values of the under- 
lying random walk (breaking ties by using the a priori order on T„). It is now quite 
obvious that, for every n > 0, the (random) empirical distribution of the random 
walk values associated with the vertices in G„ has the same distribution as X^ ■ 

Given a branching random walk BRW of the type defined above, and one of its 
vertices u, we use the notation Z(u) to denote the value of the random walk asso- 
ciated with u. The following definition is adapted from [H]. Given < v < 1 and 
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m > 1, we say that a vertex u € BRW is (m, w)— good if there is a descending path 
u =: moj ui, ■ ■ ■ , Um such that Z{ui) — Z{uo) > vi for all i € JO, mj. The importance 
of this definition comes from the two following lemmas, adapted from [TT] . 

Lemma 1. (Lemma 5.2 in [HJj Let < wi < f2 < 1. If there exists a vertex 
u E BRW at depth n such that Z{u) > V2n, then the path from the root to u must 
contain at least ^jE^^ ~ 1/(1 ~ ^i) vertices that are {m,vi)-good. 

Lemma 2. (Proposition 2.6 in ) There is a constant ijj > such that, given a 
branching random walk BRW, the probability that the root is [m, v{p) — m~'^^^)—good 
is less than exp{—ipm^^^) . 

Since Lemma [1] admits so short a proof, we reproduce it below for the sake of 
completeness. On the other hand, to give a very rough idea where the exp(—^pm^^^) 
in Lemma [2] comes from, let us just mention that it corresponds to the probability 
that a random walk remains confined in a tube of size m^/^ around its mean, for m 
time steps. Dividing the m steps into m^^^ intervals of size m^/'^, we see that this 
amounts to asking for the realization of m^^^ independent events, each of which has 
a probability of order a constant, by the usual Brownian scaling. 

Proof of LemmaUi (From pjj.) Consider a vertex u as in the statement of the 
lemma. Consider the descending path root =: xq, . . . ,Xn ■= u from the root to u. 
Let then tq := 0, and, given Tj, define inductively Tj+i := inf{j > Tj + 1; Z{xj) < 
Z{xri) +'f^i(i ~ Ti) or j = Ti + m}. Now color xq, . . . , Xn-i according to the following 
rules: if Z(xt-,+i) ^ ^i^n) + vi{j - n) and Tj+i < n + 1, then x^-., . . . ,Xri+-^-i are 
colored red. Note that this yields a segment of m red vertices, and that x^ is then 
(m,vi)— good. Otherwise, x^, ■ ■ ■ jXr^^-^-i are colored blue. Let Vred (resp. Vbiue) 
denote the number of red (resp. blue) vertices in xq, . . . ,Xn-i. Then decompose 
Z(u) into the contributions of the red and blue vertices. On the one hand, the 
contribution of red vertices is < Vred- On the other, the contribution of blue vertices 
is < Vbiue xvi + m, where the m is added to take into account a possible last segment 
colored in blue only because it has reached depth n. Writing that n = Vred + Vhiue^ 
we deduce that V2n < Vred + fi(n — Vred) + so that Vred ^ ^il^i ^ ~ — vi). 
Then use the fact that at least Vred/n^ vertices are (m, 'L'i)-good. □ 

In [llj, Lemmas [Hand [2] are used in combination with an elaborate second moment 
argument. In the present context, the following quite simple first moment argument 
turns out to be sufficient. 

Proof of TheoremUl Consider an integer r > 1, and let m := r[log(A^)^J. Let 
n > m, and let Bn denote the number of vertices in BRWi U • • - BRW at that are 
(m, v{p) — 2?7i~^/^)— good (each with respect to the BRW it belongs to) and belong 
to Go U • • • U Gn- From Lemma[T]and the definition of {Gi)i>Q, we see that the fact 
that at least one vertex in Gn has a value larger than {v{p) — m~'^/^)n implies that, 
for large n (depending on m), B^ > nm~^^^. On the other hand, Bn can be written 
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as 

(5) 

Bn := ^ 1( u is (m, v{p) - 2m-'^/^)-good)l{u G Go U • • • U G„). 

«eBRWiU-BRWiv 

Now observe that, by definition, for a vertex u of depth i, by definition, the event 
n € Go U • • • U Gn is measurable with respect to the random walk steps performed 
up to depth £, while the event that u is {m,v{p) — 2m^^/^)— good is measurable 
with respect to the random walk steps performed starting from depth > i. As a 
consequence, these two events are independent. Since the total number of vertices 
in Go U • • • U Gn is equal to A^(n + 1), we deduce from Lemma [2] and ([5]) that 
E {Bn) < N{n + 1) exp(— -0771"'^/'^). Using Markov's inequality, and letting n go to 
infinity, we deduce that 

limsupP(S„ > nm"^/^) < iVm^/^ exp(-V'm^/^). 

n— >+oo 

Now, remembering that P{Bn) = P(max using the fact 

that maxX^ < n with probability one, we finally deduce that 

limsupn~^i?(maxX^) < v{p) - m"'^^^ + Nm^^^ exp^-il^m^^^). 

n— »+oo 

Choosing r large enough in the definition of m makes the third term in the above 
r.h.s. negligible with respect to the second term, as N goes to infinity. The conclu- 
sion follows. □ 

5. The lower bound 

The proof of the lower bound on the convergence rate of vi\i{p) to Voo{p) is in some 
sense a rigorous version of the heuristic argument of Brunet and Derrida according 
to which we should compare the behavior of the particle system with N particles, 
with a version of the infinite population limit modified by a cutoff at e = 

Indeed, given a finite positive measure u on Z, let F^''-{u) be the measure defined 
by: F^'^-{v) := 2v ★ {p6i + (1 — p)Sq). This F^^' describes the evolution of the 
distribution of particles in the infinite population limit above the threshold imposed 
by the selection step. The idea of the proof of the lower bound is to control the 
discrepancy between the finite and infinite population models above this threshold. 
One important observation is that, to prove a lower bound, one does not necessarily 
have to control the number of particles at every site, but may focus instead on sites 
where the probability of finding a particle is not too small. 

We note the following two immediate properties of F^^ -. if p < v, then < 
F^^'iu), and, if g G M+, F^''-{gu) = gF^''-{v). 

5.1. Admissible sequences of measures. Throughout this section, we consider 
e > 0, < a < v{p), ^ > 1, and m > (7 := \^^^. 

We say that a (deterministic) sequence 5q =: 1^1,...,^^ of positive measures on 
Z with finite support, is (e, a, /?)— admissible, if the following properties hold: 

(i) Vi = {2py6i for < i < g; 
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(ii) for all g + 1 < i < m, vi < F^''-(m,_i); 

(iii) for all < z < m — 1, and all x E Z such that t'j(x) > 0, t'i(x) > e; 

(iv) for all q < i < m — 1, the support of I'i is contained in the interval [{v{p) — 
a){i + 1), +oo[; 

(v) Um{I^)>f3 + l. 

Note that the definition of q makes property (iv) automatic for i = q. 

Let B := {min(Xj^) < {v{p) — a)i for all 1 < i < m}. The interest of admissi- 
ble sequences of measures lies in the possibility of bounding ¥(B) from above, as 
explained in the following lemma. 

Lemma 3. Consider an {e, a, (3) — admissible sequence fo,...,fm- Let K := 
SS)^ #SMpp(i/j), and 5 := 1 — exp ^— . Then the following inequality holds: 

¥{B) < 2Kexp {-N(3-^ep6^) . 

Before proving the above lemma, we recall the following classical estimate for 
binomial random variables (see e.g. [9]): 

Lemma 4. Let n > 1 and < r < 1, and let Z be a binomial{n,r) random variable. 
Then, for allO < 5 < 1, the probability that Z < [1 — 6)nr is less than exp(— ^nr(5^). 

Proof of LemmalM For k e [[0,ml, let Ak := ClneliM^^^ - ~ '^)"'^«}- ^o^e 
that ^0 is the certain event. For < fc < m — 1, and x € Z, define Nl{x) (resp. 
Nj^{x)) to be the number of particles that are created from a particle at position x in 
during the branching step applied to , and that have a position equal to x + 1 
(resp. x). By definition, conditional upon X^, . . . , , Nlix) (resp. Nj^{x)) follows 
a binomial distribution with parameters {2X^ {x)^p) (resp. {2X^ —p))- For 
kelO,m- IJ, let Cfc := {iVf(x) > (1 - 5)2X^{x)p; x e supp(i/fc), ^ = 0, 1}. Note 
that, by definition of 5, (1 — 5)^ > /3^^ for all < < m. In view of condition (iii), 
we deduce that, for k G on A^-i, X^_^{x) > NP^^e for all x G supp(z^fc-i)- 

As a consequence. Lemma H] yields the fact that 

(6) F{Ak-i n Cfc^_i) < 2#supp(zvfc_i) exp {-N(3-'ep6^) , 

where we have used the union bound and the fact that p < 1 — p to combine the 
bounds given by Lemma H] for all the Nl_^{x), with i € {0, 1}, and x € supp(z^fc-i)- 
Now consider k G [1,^]]. On B, all the particles counted by A^^_^(A; — 1) must be 
kept after the selection step leading to X^ , since k > {v{p) — a)k. As a consequence, 
X^{k) > Nl^^{k - 1), and we deduce that 

(7) Ak-i n Ck-i nB cAk. 

Assume now that A; G + 1, mj. If B holds, we see that, according to (iv), for 
all x in the support of Uk-i, the particles counted by Nl_^{x) and N^_i{x) are all 
kept after the selection step leading to X^ . As a consequence, on Ck-i, X^ > 
(1 - (5)F^-(Xf_,l(supp(z.fc_i))), so that, on S n A^-i D Ck-i, > (1 - S)''uk, 
since X^_-^^ > (1 — S)''^^i^k-i and < F^'^'{uk-i) by assumption (ii). We deduce 
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that 

(8) n C7fc_i nBcAk. 

Now observe that, on A^, one must have X^{Z) > N{1 - 5)'^Vm{Z) > 
NP~^Um{Z) > N, a contradiction, so that Am = 0- From ([7]), ([8]), we deduce 
that F{B) < ^^Jo^ P(Afe_i n C^_i), and, using we deduce the result. □ 

Let us now relate the above results with estimates on VNip)- Define the random 
variable L := inf{l < i < m; m.m{X-^) > {v{p) — a)i}, with the convention that 
inf := m. 

Proposition 6. For all < p < 1/2, for all N > 1, 

vn{p) > {v{p) - a)(l - mF{B)). 

Proof. We define a modified branching-selection process (yJ^)„>o, composed of a 
succession of runs. Start with Lq := and Hq := 0, and i := 0, and do the 
following. 

1) Set := N5h, and k := 0. 

2) Do the following: 

{ let A; := /c + 1 and generate from the distribution pn{Z^_^, ■). } 

3) Return to 2) until k = m oi minZ^ > {v{p) — a)k + Hi. 

4) Set Lj_|_i := Li + k and -ffi+i := min Z^ 

5) Set (yi^,...,yiy^^_,):=(Zo^...,zf_J 

6) Let i := i + 1 and return to 1) for the next run. 

One may describe the above process as follows: starting from a reference posi- 
tion Hi and a reference time index L,, a run behaves like the original branching- 
selection process until either m steps have been performed or the minimum value 
in the population of particles exceeds the reference position by an amount of 
at least {v{p) — a) times the number of steps performed since the beginning of 
the run. Then the current population of A*" particles is collapsed onto its min- 
imum position, the reference position is updated to this minimum position, and 
the time index to the current time, and a new run is started. Our modified 
process has a natural regeneration structure yielding the fact that the sequences 
(Lj+i — Lj)j>o and (-ffi+i — Hi) are i.i.d. The common distribution of the -Lj+i — L, 
is that of L, while the common distribution of the -ffi+i — Hi is that of minX^ . 
From the fact that minY^^ > Hi, it is easy to deduce that, with probability 

one, lim„^+oo 'T'"^ min = and, since < n~^Yj^ < 1, we also 

have that lim„^+oo n"^E(minl^) = Now, by definition, minX^ > 

{v{p) - a)Ll(S=), so that E(minXf ) > {v{p) - a)(E{L) - E{L1{B))). Using the 
fact that 1 < L < m, we obtain that ^^'^^^^ > {v{p) — «)(! — mF{B)). 



Now, it should be intuitively obvious that the modified process )„>o is in some 
sense a lower bound for the original process {X^)n>o, since we modify the original 
dynamics in a way that can only lower positions of particles. It is in fact an easy 
consequence of Proposition [1] that, for all n, the distribution of Y^^ is stochastically 



10 



JEAN BERARD 



dominated by that of Xj^ . A consequence is that E(miny^) < E(minX^).The 
result follows. □ 

5.2. Construction of an admissible sequence of measures. Let A denote an 
integer > 4, and m > q. Then let am ■= [m^^^l, Cm ■= [rn^^^l, Sm ■= L 2(i-"(p)) J- 
Define dm by dm (fc) := for /c E [1, Sml, and (irrt(fc) ■= v{p)k+am ^or k £ w.]]. 
Define gm by gm{k) ■= kfor k € [1, (k) ■.= v{p){k+l) fovke ISm + l 

and gm{k) ■= v{p)k — Aam for k € [m — Cm + 1, m^- 

Then define a sequence of measures (t'j)ig|o,m] ^ follows. Let iSi)j_^^Q„^j 
denote a simple random walk on Z starting at zero, with step distribution p5i + (1 — 
p)6o, governed by a probability measure P, then let 

Ui{x) := 2'P [gm{k) < Sk < dm{k) for all k G ^0, ii,Si = x\. 

The main result in this section is the following: 

Proposition 7. For large enough A, there exists x(^) > such that, for all large 
enough m, the above sequence is {exp{—xiA)m^^^ ,2Am''^^^ ,2008)-admissible. 

We need to establish several results before we can prove the above proposition. 
First, consider the modified probability measure P defined by 

dp [v{p)\^^ (l-v{py'-^^ 



dP \ P J V 1 ~ P 

With respect to P, (»S'j)jg|o^m] is a simple random walk on Z starting at zero, with 
step distribution v{p)Si + (1 — v{p))5o. 

We now rewrite ^'i(x) in terms of this change of measure. To this end, introduce 
the compensated random walk defined by Si := Si — v{p)i, let also gm{k) := gm{k) — 
v{p)k and dm{k) := dm{k) - v{p)k. Finally, let 7 := ^(p)/(jl^(p)) , and note that 
7 < 1 since p < v{p). After a little algebra involving the definition of v{p) in terms 
of A, we obtain the following expression: 

(9) Ui{x) = E (^gm{k) <Sk< dm{k) for ah k G JO, ij, Si 

where E denotes expectation with respect to P. From the definition of gm and 
dm, we see that, for large m, the only values of Si that contribute in the above 
expectation are < m^^^. As a consequence, 

(10) Ui{x) > 7"'^'P [gm{k) < Sk < dm{k) for ah k € lO, il S^ = x 
Lemma 5. For some constant C,i > 0, as m goes to infinity, 



P 



9m{k) < Sk < dm{k) for all k G [[0,m| > exp(-Cim^/^ 



We need an elementary lemma before the proof of Lemma [5j 

Lemma 6. Consider a random walk {Zi)i>Q on M, defined hy Zi := Zq + £i + ■ ■ ■ +£i 
for i > 1. Assume that the increments £i are i.i.d. with respect to some probability 
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measure Q and satisfy E{ei) = and < Var{ei) < +oo. Then there exists A > 
such that, for all m large enough, on the event that a-m/'i ^ Zq < 2am /3, 

(11) Q[v < Zi<am; « G (0, Cmj, am/3 < Zc^ < 2am/3| Zq] > A, 

(12) Q [am/4 <Zi< 3am/4; i G [[0,Cm]]| ^o] > A. 

Proof of Lemma 0. We use the convergence of the distribution of the random process 
(A^r)iG[o,i] defined by := 0, M™^^ = a-\ei + ■ ■ ■ + e^) ior i e [[l,Cml, and 
interpolated hnearly on each [i/cm,{i + l)/cm] towards the Brownian motion (on 
the space of real- valued continuous functions on [0, 1] equipped with the sup norm). 
An easy consequence is that there exists Ai > such that, for all large m, 

Q [-am/7 < Zi- Zo < am/7; i G [1, Cml, < Zc^ - Zq < am/7] > Ai. 

This estimate proves (jlip when Zq belongs to [am/3, am/2]. A symmetric argument 
works when Zq belongs to [am/2, 2am/3]. The proof of (fT2]) is quite similar. □ 

Proof of LemmalM First note that, for all large enough m, am/3 < Sm — v{p)sm < 
2am/3. Then, 

P [gmik) <Sk< dm{k); k G [0,sml] = v{py^. 

Divide the interval [[sm + into consecutive intervals Ij, j G [[1,/im]]) 

where each Ij for j G [[1,/im — l] is of the form Ij := [[sm + 1 + (i — 
l)cm,Sm + 1 + jcmj, while the last interval is := [[sm + 1 + (hm - l)cm,'m-]], 
whose length is less than or equal to Cm- For i G Jl,/im]], let bm,j-i and 
bm,j be defined by Ij = [[6mj-i, &mj]]- Now, for i G [[1,/im - 11, define 
the event F, := < Sk < Om] k G Ij, am/3 < Zb^ ^ < 2am/3|, and let F^^ := 

{am/4 <Sk< ?>am/^] A: G 

It is easily checked that, given that Ss^ = Sm, 

f] TiC fl {gm{k) <Sk< dm{k)}. 

je[i,/iml fce[Sm+l,m] 

Thanks to Lemma [6] and to the Markov property of Z with respect to P, we deduce 
that 



P 



9m{k) <Sk< Mk); k G [0,mj > vipY'- X^" 



(We use exactly Lemma [6] for intervals /j with i G p,/im — ll, while, for i = km, 
we use that fact that the length of Ij is < Cm, whence the fact that the conditional 
probability of F^^ given Sq, . . . , S^^ h,n-\ larger than or equal to the probability 
appearing in (jl2p in Lemma [6l) Using the fact that hm ~ m^/^ for large m, the 
conclusion follows. □ 

Lemma 7. There exists C2(^) > such that, as m goes to infinity, 

m({i^i{xy, f G [0,m - ll, Ui{x) > 0} > exp{-C2{A)m^^^). 

We need the following lemma before giving the proof. 
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Lemma 8. Let p, cr € M, with p + 1 < a, v g]0, 1[ and, let £ be an integer such that 
a+£v < p+£ and p+vi > a. Then, for all x G Zn[p, a], and all y € Zn[p+w£, a+vl], 
there exists a sequence x =: xo,xi, . . . ^xi := y such that Xj+i — Xj € {0, 1} for all 
i e [0, ^ - 1|, and p + vi < Xi < a + vi for all i G [0, Ij. 

Proof of LemmalM Consider x £ Zr\[p,a], and y G Z n [p + vi, a + vi]. Define 
inductively the sequence (rj, /ii)i>o as follows. Let tq := x, ho := x. Our assumption 
that p + 1 < a guarantees that x or x + 1 belongs to [p + v, a + v]. If x + 1 belongs 
to [p + v,a + v], then let d := 0. Otherwise, let d := 1. Then, consider i > 1. If 
i + d is even, let Tj := max{j G Itj-i + 1, +oo[[; hi-i > p + vj}, and let hi := 
If i + d is odd, let Tj := max{j G lu^i + 1, +oo[[; + j — rj_i < cr + f j}; and let 
hi := /ij-i + Tj — Ti^i. 

The fact that and p + 1 < a and f g]0, 1[ guarantees that every term in the 
sequence is finite. Define the sequence {zk)k>o by -Zfe := K-i for k G [[Tj_i,rj]] when 
i + d is even, and := /ii-i + A; — rj_i for k G [[rj_i,rj]] when i + d is odd. Our 
assumptions yield the fact that p + vk < < a + vk for all > 0. Now note that, 
il y = Zmi the path zq,. . . ,Zm solves our problem. If y > Zm, the assumption that 
p + £ > CT + v£ plus elementary geometric considerations show that there must be 
some A; G [0, m — 1| such that y — x^ = i — k. From the path xq, . . . ,Xk, add +1 steps 
and stop at length i. The corresponding path solves the problem. Now, if y < Zm, 
our assumption that a < p + v£ plus elementary geometric considerations show that 
there must he a k £ [[tk_i, m — 1| such that x^ = y. From the path xq, . . . , Xfc, add 
+0 steps and stop at length £. The corresponding path solves the problem. □ 

Proof of Lemma^ In the sequel, let x^ := k for k G JO, s^J, and note that, for large 
enough m, gmik) < Xk < dm{k) for all k G [0, s^J. 

Let a > max{v{p)~^ , (1 - v{p))~^), and let fm ■= [aam]- For i G [s-m + l,Sm + 
2/m]], and x in the support of fj, there must by definition be a sequence Sm =■ 
Xs^, ■ ■ ■ , Xi := X such that Xj+i — Xj G {0,1} and gm{j) < < dm{j) for all 
j £ [[•^rrn'^l- As a consequence, for every i G [[0, + 2/m]] and x G supp(fj), 

(13) P[5fc = xfc; kelO,il Si = x]>{min{v{p),l-v{p)y^+^^^. 

Now assume that i G [[sm+2/m+l, rn,— CmJ, and consider the distribution of Si-j^ 
with respect to P, conditional upon gm{k) < Sk < dm{k) for all A: G [0,? — /mj. 
This probability distribution is concentrated on the set Z n [gmii — fm),dm{i — fm)], 
which contains at most elements. Consequently, there exists u G Z n [gmii — 
fm),dm{i - fm)] such that 

P[Si 

—fm ~ ^1 9m (k) < Sk < dm{k); k £lO,i- fmj] > l/a-m. 

Now let X belong to the support of Vi. By definition, x G Z n [5m(^)) c?m(^)]- In view 
of the definition of a, we see that, for all m large enough, we can apply Lemma [8] to 
obtain the existence of a sequence u =: Xi-f^, . . . , Xj := x such that Xj+i—Xj G {0, 1} 
and gm{j) < Xj < dm{j) for all j G p — fm,^- Thanks to the Markov property of 
{Sk)k>o with respect to P, we deduce that 

P [gmik) <Sk< dm{k); k G [[0,i]], Si = x] 
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is larger than or equal to 

P [9m{k) < Sk < dmik), kelO,i- /„]]] (l/a^)(min(t>(p), 1 - 
From Lemma [5l we have that, for m large enough, 

P[9m{k) <Sk< dm{k), k e [[0,ml] > expi-Cim^/'^), 
so that, for every i € [[sm + 2/m + 1, m — Cmj, 

(14) P [graik) <Sk< dra{k), k G [0, ij, = x] > exp(-Cimi/3) 

x(l/am)(min(t;(p),l - vip))^"^ . 

For i € Im — Cm + ^, rn — Cm + 2(^ + I)/™]], any x in the support of i^i(x) is such 
that there exists u in the support of Vm-cn^ and a sequence u =: Xm-c„^, ■ ■ ■ ,Xi := x 
such that Xj+i — xj £ {0, 1} and QmU) ^ ^ dm{j) for all j G Jm — Cm,i]]- As a 
consequence, 

P bm(fc) < 5fc < dmik); k G [[0,^1, 5, = x] 
is larger than or equal to 

P [gmik) < Sk < dmik); k elO,m- Cmh Sm-c^ = u] 

x(min(t)(p), 1 - 

Using (dH), we deduce that, for every i G — Cm + ^, m — Cm + 2(^4 + l)/m]]) 

(15) P [gmik) <Sk< dmik); k G [[0,i]], 5, = x] > exp(-Cimi/3) 

X (l/a^)(min(t;(p), 1 - ^;(p))(2(^+i)+i)/- . 

Then an argument quite similar to that leading to (|14p yields that, for any i G 
lm-Cm + 2iA + l)fm + hml 

(16) P b™(A:) < 5fc < dmik); k G [[0,^1, 5, = x] > expi-dm^/'^) 

x(l/((A + l)a^ + l))(min(t>(p), 1 - vip))^^+^^f"^ . 
The conclusion follows by using ((l0|) and ([13]), ([II]), ([IS]), ([H]). □ 
Lemma 9. There exists (piA) > such that, as m goes to infinity, 

P [gmik) <Sk< dmik); k G I0,m]], Sm < -(^/2)a^] > (/)(A) exp(-Cimi/3). 

We shall need the following elementary lemma. 

Lemma 10. Consider a random walk (Zj)j>o on M, defined by Zi := Zo+ei + - • -+£1 
for i > 1. Assume that the increments £i are i.i.d. with respect to some probability 
measure Q and satisfy -E'(ei) = and < Variei) < +00. Then there exists 
4>iA) > such that, for all m large enough, on the event am/4 < Zq < 2>am/^, 

Q[-Aam <Zi<am] ie U,Cml Zc^ < -iA/2)am\ Zq] > (piA). 

Proof of Lemma \10l We re-use the notations of the proof of Lemma [H The only 
point is to note that, by convergence to the Brownian motion, there exists (piA) > 
such that for all large m, 

Q [-Mam/4.<Zi - Zo<am/4.; i G [[l,c™]], Z^^ - Zo<-(A/2 + 8/4)0^] > 0(A). 
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The result follows easily (using the fact that A is assumed to be > 4). □ 
Proof of Lemma O From the proof of Lemma [SJ we see that 



P 



gm{k)<Sk<dm{k);k e [0, m - Cml, am/4<5m_c„<3am/4 > exp(-Cim^/^) 



Then, by Lemma [TOl as m goes to infinity, 



P 



9m{k)<Sk<dmik);k £ Im - Cm + l,rr4, Sm<-Aam/2 am/4<S'm-c„<3am/4 



is larger than or equal to (j){A). The conclusion follows. □ 

Proof of Proposition Properties (i) and (ii) are immediate consequences of the 
definition. Property (iii) is a direct consequence of LemmaO letting x(^) •= C2(^)- 
From the definition, (k) > v{p){k + 1) for all k^lq , ni — CmJ- Then, given A, for 
all large enough m, it is easily checked that gm{k) > {v{p) — 2Am~'^^^){k + 1) for all 
^ — H~ l^^l" This yields Property (iv). As for Property (v), consider (|9|) . 
Clearly, since 7 < 1, 

i^m{^) > 7"(^/2)»-P ^g^{k) <Sk< dm{k); k G [0,m]], 5^ < -{A/2)a, 

From Lemma [U the probability in the r.h.s. of the above expression is > 
(p{A) exp(—(^im^^^), so that, choosing A large enough, the term ^~(^/2)am f^omi- 
nates for large m. As a consequence, for such an A, Umi'^) > 2008 + 1 as soon as m 
is large enough. □ 

Proof of Theorem\M Consider a parameter 9 > 0, and let m depend on N in the 
following way: m := [01og(A^)J^. For N large enough, we can apply Proposition [71 
and use Lemma [3l which states that 

r{B) < 2Kexp {-N(3-^ep6'^) . 

It is easily checked from the definition that K < {(A + l)ar« + l)™-. Note also that 
5 ~ log(/3)/m for large N. Finally, e = exp(— x(^)"i^^^). Choosing 6 small enough, 
we check that mF{B) « m~^/^ as A'^ goes to infinity. Proposition [6] then implies 
that VNip) > v{p) - 2Am"2/3(l + o(l)). The result follows. □ 

6. The case 1/2 < p < 1 

In the case l/2<p<l, it turns out that Voo{p) = v{p) = 1, which makes it much 
easier to obtain estimates. 

6.1. Upper bound when p = 1/2. It is easily checked that, for all m > 0, the 
number of particles at position exactly m after m steps, that is, X^{m), is stochas- 
tically dominated by the total population at the m— th generation of a family of 
N independent Galton- Watson trees, with offspring distribution binomial(2, 1/2). 
This corresponds to the critical case of Galton- Watson trees, and the probability 
that such a tree survives up to the m— th generation is < cm~^ for some constant 
c > and all large m. As a consequence, for large enough m, ¥{X^{m) > 1) < 
E{X^{m)) < cNm-K 
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On the other hand, we have by definition that m ^Km.ax{X^) < 1 — 
^P(X^(m) = 0). Choosing m := AN, where ^ > 1 is an integer, we see that, 
for large N, m-^Emax(X^) < 1 - 1/AN{1 - c/A). The upper bound in 1^ follows 
by choosing A > c. 

6.2. Lower bound when p = 1/2. Given m > 1, define U := inf{n G 
[[l,m]]; (n) < 2N/3}, with the convention that inf0 := m. Let D denote the 
event that minX^ < U — 1. 

Using an argument similar to the proof of Proposition [6l with U and D in place 
of L and B respectively, we deduce that 

(17) VN{p)>l-^^-mF{D). 

The lower bound in ([3]) is then a consequence of the two following claims. 

Our first claim is that there exists c > such that ¥{D) < exp{—cN) for all 
large A^. Recall the definition of N^{x) from the proof of Lemma [3l and choose 6 
small enough so that (1 - 5)4iV/3 > A^. It is easily seen that D C {NIj_^{U - 1) < 
(1 - 6)2N/3} U {N^_i{U - 1) < (1 - 6)2N/3}. Now, by definition, one has that 
^U-iiU — 1) > 2N/3, so that, by the bound of LemmalU conditional on X^_^, the 
probabilities of N}j_-^{U - 1) < (1 - 5)2N/3 and N^_i{U - 1) < (1 - 5)2N/3 are 
both < exp(-c((5)7V). The bound on P(L») follows. 

Our second claim is that, for small enough e > 0, with m := [eA^J, there ex- 
ists c(e) > such that E(D) > c{e)N for all large A^. To prove it, introduce the 
Markov chains (Vfc)A;>o and {Zk)k>o on N, defined as follows. First, Vq G N, and, 
given Vq, . . . , Vfc, the next term V^+i is the minimum of A^ and of a random vari- 
able with a binomial(2T4, 1/2) distribution. On the other hand, Zq E N, and, given 
Zq, . . . , Zk,the distribution of Zk+i is binomial(2Zfc, 1/2). Observe that the sequence 
{X^ {n))n>o is a version of V started at Vq := A^. Now, it is easily seen that, given 
two starting points x, y € N such that x < y, one can couple two versions of V start- 
ing from X and y respectively, in such a way that the version starting from y is always 
above the version starting from x. As a consequence, U stochastically dominates the 
random variable T := inf{n € [1, mj; Vn < 2A^/3} (again with inf := m), where V 
is started at Vq := [3A'^/4J. Then observe that the distributions of V and Z started 
with Vq := Zq := [3A^/4J , considered up to the hitting time of [A^, +oo[[, coincide. As 
a consequence, the probabilities of the events {sup^g|Q ,„j |Vfc— [3A^/4J | > A^/16} and 
{supj;,g|o^m.] — [3A^/4J I > A^/16} coincide, and the first of these two events implies 
that T = m. Now, (.^fc)fc>o is a martingale, so that, by Doob's maximal inequal- 
ity, P (supfcg|o,H \Zk - [3N/4\\ > N/w) < E{Z^ - \3N / A\f {N / Then, it 

is easily checked from the definition that E{Z'^_^-^\Zk) = Z'^ + Zk/2, and, using again 
the fact that {Zk)k>o is a martingale, we deduce that E{Zm - l3N/4\f < mN/2. 
As a consequence, we see that, choosing e > small enough, we can ensure that 

P (Efcelo,Le7VJ] \Zk - [3N/4:\ \ > N/loj < 1/2008 for all large A^. For such an e, and 
all N large enough, we thus have that F{U = m) > P{T = m) > 1/2008. The 
conclusion follows. 
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6.3. Upper and lower bound when 1/2 < p < 1. As for the upper bound, 
observe that asking all the 2N particles generated during the branching step to 
remain at the position from which they are originated has a probability equal to at 
most (1 — p)^^ , so that E(maxX„) < n(l — (1 —p)^^)- As for the lower bound, 
observe that, starting from N particles at a site, the number of particles generated 
from these during a branching step and that perform +1 random walk steps has 
a binomial(2A^,p) distribution, whose expectation is 2pN , with 2p > 1. Using the 
bound stated in LemmalU we see that the probability for this number to be less than 
is < exp(— cA^) for some c > 0. As a consequence, IE(minX„) > n(l — exp(— cA^)). 

7. Concluding remarks 

Remark 1. Can we derive a reasonably simple explanation of how the log(A^)^^ 
arises, based on the mathematical proofs presented above ? Broadly speaking, the 
key point in both the upper and the lower bound seems to be the following (we re-use 
some notations from Section\^^: consider a large integer m and look for a scale A 
such that 

(18) 2'"P(|5i|,...,|5„|aA)xl/Ar. 

In view of the change of measure, and of the fact that, by Brownian scaling, 
P{\Si\, . . . ,\Sm\ C)C A) X exp(— m/A^), we see that there are two factors involved 
in the above probability: 7^'^, and exp(— m/A^). Equating the exponential scales of 
these two factors yields A oc m^^^ , and il8\) then implies that m oc log(A^)^, whence 
an average velocity shift over the m steps of order A/m oc log(A^)^^. 

Remark 2. What we have proved is that the order of magnitude of vn{p) — Vooip) is 
indeed log(A^)^^. It would of course be quite interesting to get more precise asymp- 
totic results for this quantity, since there is at least compelling numerical evidence 
that vn{p) — Voo{p) ~ c{p) log(A^)~^ for some constant c. 

Remark 3. Both the upper and lower bound presented here relic upon controlling 
the behavior of the particle system for time intervals of length m oc \og{Nf. This 
is the same order of magnitude as the one observed for the coalescence times of 
the genealogical process underlying the branching- selection algorithm, from empirical 
studies and heuristic arguments (see e.g. [3]J. Although we do not know how to 
establish a rigorous relationship between these facts, this at least provides another 
indication that the log(A^)^ time scale is particularly relevant for the study of these 
systems. 
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